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The Nambu - Jona-Lasinio model in 2+1 dimensions is simulated for non-zero baryon 
chemical potential with a diquark source term. No evidence for a BCS condensate or 
gap is seen at high density; rather, critical behaviour with novel exponents is observed, 
suggesting that 2d superfluidity as first described by Kosterlitz and Thouless is realised, 
but with the universality class determined by the presence of relativistic fermions. 

1. INTRODUCTION 




Figure 1. The BCS instability in degenerate matter. 



There has been intense recent speculation on the form of the QCD phase diagram 
at small temperature T but large baryon chemical potential /i. The conventional view 
is that for sufficiently large yU, chiral symmetry is spontaneously restored, resulting in 
a degenerate system of relativistic quarks known as quark matter, with Fermi energy 
Ep = /i/3. Since the qq interaction can be attractive due to gluon exchange and/or 
instanton effects 0, however, this simple picture may have a BCS instability with respect 
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to condensation of diquark pairs at the Fermi surface, resulting in a color superconducting 
ground state [Q, and the formation of an energy gap A 7^ to the first excited state of 
the system (see Fig. |ID. Model calculations suggest that A could be as large as 100 
MeV, comparable with the constituent quark scale. 

It is clearly desirable to examine the various scenarios in non-perturbative lattice QCD 
simulations, where systematic control is at least in principle possible. Unfortunately, this 
has so far proved impracticable because the Euclidean functional measure is no longer 
positive definite for 7^ 0, rendering the usual importance sampling methods ineffective. 
There are, however, simpler four-fermion models where Monte Carlo methods can be 
applied and yield a plausible description of degenerate matter p. In this talk I describe 
numerical studies of the Nambu - Jona-Lasinio (NJL) model in 2+1 spacetime dimensions; 
the Lagrangian in continuum notation reads 

C = i,{^ + /i7o + m)V^ - |- [(7/^^)2 - [i^^^fijf] . (1) 

The model in 3+l(i has a long history as an effective description of strong interaction 
physics 0. For coupling stronger than some critical g1 the SU(2)i®SU(2)/j®U(l)B 
global symmetry present for quark mass m — » spontaneously breaks to SU(2)/®U(l)s, 
accompanied by the generation of a constituent quark mass S = g'^i^ipip). The behaviour in 
2+ld has the same pattern, except that this time there is an interacting continuum limit 
dX g"^ gli S/A — >■ Once a chemical potential is introduced, for low T an expansion 
in the inverse number of flavors 1/A^/ predicts to remain roughly unchanged out to a 
critical /ic — S(O), whereupon chiral symmetry is abruptly restored in a strong first-order 
transition [0. At the same point the baryon density ub = (ip'joip) rises from zero to a 
non-zero oc fi'^6{fi — fic), consistent with filling a two-dimensional Fermi ball of radius 
jj,. This is confirmed by lattice simulations 0; in particular, unlike other simulations with 
a real measure, there is a clear separation of scales between fi^ and m-„ [Q . 

2. THE DIQUARK SECTOR 

The question to consider is whether for fi > fic the baryon number U{1)b symmetry 
is spontaneously broken by a diquark condensate (qq) 7^ 0. Since the NJL model is not 
a gauge theory, this leads not to superconductivity but to the associated phenomenon of 
superfluidity. In fact, numerical studies ||^ reveal that the preferred channel for pairing 
in this regime is (in the notation of staggered fermions) via the scalar SU(2)2,®SU(2)r 
singlet x*'^T2X. To find unambiguous evidence for BCS condensation on a finite system. 



however, the correct procedure ||T0[ is to add to £ a diquark source term 

J±qq± = J±{x''r2X ± Xrax*")- (2) 

It is then possible to measure the diquark condensate using methods similar to those used 
in lattice studies of chiral symmetry breaking: 

, V IdlnZ 

= (3) 

together with the associated susceptibilities 

X± = I](9^±(0)^9±(^))- (4) 
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In a conventional BCS condensation the + describes a "Higgs" mode while the - is a 
"Goldstone" , which is constrained by a Ward identity: 



x-iJ-= 0) 



(5) 
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• 1^=0.2 
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Figure 2. ln(gg+) vs. Inj. 



Figure 3. vs. j. 



In collaboration with Biagio Lucini and Susan Morrison I have simulated a lattice NJL 
model including the source term (|]) [|lT|. Our parameter choice yields a /ic — 0.65. Fig. |^ 
plots extrapolated to the zero temperature — > cxd limit; for \i < \Xf, the data 

support a linear relation between and j, consistent with a U(l)B-symmetric ground 

state. For /i > by contrast, the data suggest 

oc j" (6) 

with a = a{fi) falling in the range 0.2 - 0.3 for the fi values studied. Evidence for power 
law behaviour is reinforced by considering the susceptibility ratio \x+/X-\- eqns. (|^) and 
(|^) together imply that this ratio should take the constant value a in the high density 
phase, and the plateaux of Fig. |^ for /z = 0.8, 0.9 support this. 

The simulation data suggest that the high density phase > /ic is critical, characterised 
by continuously varying exponents = and //(/i) defined via 

(gg(O)gg(x)) oc xe R^ (7) 

and is thus qualitatively similar to the low-T phase of the 2d XY model. If we write 
qq{x) = (t)oe^^'^^\ then long range order is washed out by IR fluctuations of 6, but long- 
range phase coherence persists via (J^). This is precisely what gives rise to persistent 
currents and hence superfluidity in 2d systems. The supercurrent Jg is related to qq via 

Js = KsVe. (8) 
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The only way to change the circulation n = § Js.dl around a periodic volume is to create 
a vortex - anti-vortex in the {6} configuration and translate one of the pair around the 
universe in the perpendicular direction until they reannihilate, in so doing incrementing 
the quantised k by 2ttKs/L. The energy needed increases logarithmically with L, implying 
that the circulation is metastable 
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3. THE SPIN-i SECTOR 
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Figure 4. Amplitudes for hole {A), Figure 5. Dispersion relation E{k) for both 

particle {B) and anomalous (C) propagation, free and interacting fermions. 



The conclusions of the previous section beg the question of why the observed 5^3 
- 5 is not consistent with the 2d XY value S{T) > 15, in contrast to a recent study 
of superfiuidity in the attractive Hubbard model [|1^]. Further insight is gained from 
studying the spin-^ sector via the fermion propagator, which for j 7^ contains both 
normal {q{0)q{x)) and anomalous {q{0)q{x)) components. To probe the Fermi surface, 
data at non-zero momentum k are needed. In the normal sector a sharp transition between 
forwards Ae~^^ and backwards 5e~'^^^*~*^ propagation at a rather well-defined momentum 
/c/tt ^ 0.28 on a 32'^ system at yU = 0.8, as shown in Fig. ^ We interpret this as a crossover 
between hole and particle excitations at the Fermi surface. The anomalous propagators 
signifying particle-hole mixing peak at the same momentum, but appear to vanish as 
j^O. 

The results for the mass gap are shown in Fig. with hole energies plotted as negative. 
The resulting curve is the quasiparticle dispersion relation; its detailed form indicates a 
Fermi momentum kp slightly less than fi and a Fermi velocity Pf = dE/dk\k=kp — 0.7, 
somewhat less than the speed of light. Both results differ from the corresponding free- 
field values, but are consistent with a relativistic normal Fermi liquid with repulsive 



interaction between quasiparticles with parallel momenta [jl^. Most importantly, there 
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is no signal for a gap A 7^ 0. The origin of the non-standard critical behaviour may 
therefore be attributed to massless fermions; the presence of a Fermi surface leads to a 
new 2d universality class. This situation should be contrasted with dimensional reduction 
in systems with T > ; because the fermi degrees of freedom in this case acquire masses 
ttT and decouple, the universality class is that of the spin model with the corresponding 
global symmetry. 



4. CONCLUSION 

Both [qq) and A vanish in the limit j — > implying that the conventional BCS descrip- 
tion appropriate for, say, superfiuid ^He, does not apply in this case. The most economic 
hypothesis explaining the results of Figs. 2 - 5 is that the NJL2+1 model in its high density 
phase describes a gapless relativistic system with superfluidity realised in the way first 
suggested by Kosterlitz and Thouless for thin films of ^He. Further numerical confirma- 
tion for this picture would come from observation of current quantisation in the presence 
of a source j{x) with built-in winding number. The more urgent problem, however, is 
to extend the calculations to NJL3+1 and test the exciting scenarios predicted in 
beyond the self-consistent approximation. 
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